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Abstract 

Background: Although some studies have indicated that climate changes can affect Pinus koraiensis mixed forest, the 
responses of composition and structure of Pinus l<oraiensis mixed forests to climatic changes are unknown and the key 
climatic factors controlling the composition and structure of Pinus l<oraiensis mixed forest are uncertain. 

Methodology/principal findings: Field survey was conducted in the natural Pinus l<oraiensis mixed forests along a 
latitudinal gradient and an elevational gradient in Northeast China. In order to build the mathematical models for simulating 
the relationships of compositional and structural attributes of the Pinus l<oraiensis mixed forest with climatic and non- 
climatic factors, stepwise linear regression analyses were performed, incorporating 14 dependent variables and the linear 
and quadratic components of 9 factors. All the selected new models were computed under the +2°C and +10% 
precipitation and +4°C and +10% precipitation scenarios. The Max Temperature of Warmest Month, Mean Temperature of 
Warmest Quarter and Precipitation of Wettest Month were observed to be key climatic factors controlling the stand 
densities and total basal areas of Pinus l<oraiensis mixed forest. Increased summer temperatures and precipitations strongly 
enhanced the stand densities and total basal areas of broadleaf trees but had little effect on Pinus l<oraiensis under the +2 C 
and +10% precipitation scenario and +4''C and +10% precipitation scenario. 

Conclusions/significance: These results show that the Max Temperature of Warmest Month, Mean Temperature of Warmest 
Quarter and Precipitation of Wettest Month are key climatic factors which shape the composition and structure of Pinus 
l<oraiensis mixed forest. Although the Pinus l<oraiensis would persist, the current forests dominated by Pinus l<oraiensis in the 
region would all shift and become broadleaf-dominated forests due to the dramatic increase of broadleaf trees under the 
future global warming and increased precipitation. 
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Introduction 

Pinus koraiensis is a pine species endemic to the region including 
East Russia, Korea, Japan and Northeast China [1,2], which 
provides edible seeds and very useful timber [2], and usually forms 
mixed forest with various broadleaf trees [3] . The Pinus koraiensis - 
broadleaf mixed forest is one of the major vegetation types and the 
typical conifer - broadleaf mixed forest and plays important role in 
carbon cycling in Northeast China [3,4]. 

A field study of the North East China Transect (NECT) has 
observed the shrinkage of Pinus koraiensis patches and the expansion 
of some broadleaf tree species from 1986 to 1994, and indicated 
that the reason for the shrinkage of Pinus koraiensis patches and the 
expansion of some broadleaf tree species might be the effects of 
climate change [5]. Some simulation studies also showed that 
climate changes could affect Pinus koraiensis or its mixed forest 
[6,7]. For example, by using an ecological information system 
(GREEN), Xu and Yan (2001) predicted that the potential 
distribution area of Pinus koraiensis would shift Northward under 



the HadCM2 climatic change scenario and under the average of 
another five climatic change scenarios (GISS, NCAR, OSU, 
UKMO and MPI) [6]. By linking a forest gap model (LINKAG- 
ES) with a landscape model (LANDIS), He et al. (2005) predicted 
that the broadleaf trees would overtake Pinus koraiensis to become 
broadleaf forests in Changbai Natural Reserve of Northeast China 
under the CGCM2 climate change scenario [7]. However, the key 
climatic factors controlling the composition and structure of Pinus 
koraiensis mixed forest are uncertain. 

Linear or quadratic regression models can be built to simulate 
the relationships of compositional and structural attributes of plant 
communities with climatic and non-climatic factors, and such 
studies have been conducted in grasses, tropical trees, Pinus contorta 
and Pinus sylvestris [8,9,10,11,12]. To better understand the 
compositional and structural changes of plant communities, it is 
necessary to recognize the effects of natural succession with age, 
which strongly affect the compositional and structural changes of 
plant communities even without evidence of climatic change 
[13,14,15]. Usually, the mean diameter of Pinus koraiensis reflects 
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Figure 1. Location of study sites in Northeast Cliina. 

doi:1 0.1 371/journal.pone.00971 92.g001 



Table 1. 


Location and elevation of Pinus koraiensis mixed forests in 


different sites and plots. 






Site 


Plots 


Plot Latitude (°N) 


Plot Longitude (°E) 


Plot Elevation (m) 


Kuandian 


Three 30 m * 20 m plots 


40.9123 


1 24.7883 


930 


Dongsheng 


Three 30 m * 20 m plots 


44.4113 


128.1233 


804 


Llangshui 


Three 30 m * 20 m plots 


47.1662 


1 28.8843 


445 


Shengshan 


Three 20 m * 20 m plots 


49.4780 


126.7782 


504 


Changbai A 


Three 30 m * 20 m plots 


42.3993 


1 28.0946 


755 


Changbai B 


Three 30 m * 20 m plots 


42.4001 


128.0952 


759 


Changbai C 


Three 30 m * 20 m plots 


42.3742 


128.0852 


851 


Changbai D 


Three 30 m * 20 m plots 


42.2655 


128.1476 


960 


Changbai E 


Three 30 m * 20 m plots 


42.2270 


128.1743 


1058 


Changbai F 


Three 30 m • 20 m plots 


42.1914 


128.3117 


1076 


doi:l 0.1 371 /journal.pone.00971 92.t001 
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Table 2. The 23 Variables and their abbreviations used in the study. 





Variable 


Abbreviation 


Max Temperature of Warmest Month (°C) 


TempWarmestMonth 


Min Temperature of Coldest Month {"C) 


TempColdestMonth 


Mean Temperature of Warmest Quarter ( C) 


TempWarmestQuarter 


Mean Temperature of Coldest Quarter (X) 


TempColdestQuarter 


Precipitation of Wettest Month (mm) 


PrecipWettestMonth 


Precipitation of Driest Month (mm) 


PrecipDriestMonth 


Precipitation of Wettest Quarter (mm) 


PrecipWettestQuarter 


Precipitation of Driest Quarter (mm) 


PrecipDriestQuarter 


Mean DBH of Pinus koraiensis (cm) 


MeanDBHpinus 


Total Basal Area of Pinus koraiensis [m^ ha~^} 


BasalAreapinus 


Stand Density of Pinus koraiensis (No. ha~^) 


StandDensitypinus 


Total Basal Area of /Acer and Fraxinus (m'^ ha~^) 


Basal Area/\(-er-Fraxinus 


Stand Density of >Acer and Fraxinus (No. ha~^) 


StandDensityAcer-Fraxinus 


Total Basal Area of Quercus (m^ ha~^) 


BasalAreaquercus 


Stand Density of Quercus (No. ha~^) 


StandDensityquercus 


Total Basal Area of Abies and Picea (m^ ha~^) 


BasalAreaAbies-Picea 


Stand Density of Abies and Picea (No. ha~^) 


Stand Den sityAbies-Picea 


Total Basal Area of Larix (m^ ha~^) 


BasalAreaLarix 


Stand Density of Larix (No. ha~^) 


StandDensityLarix 


Total Basal Area of Ulmus (m^ ha"^) 


BasalAreauimus 


Stand Density of Ulmus (No. ha~^) 


StandDensityuimus 


Total Basal Area of Betula (m^ ha"'') 


BasalAreaeetuia 


Stand Density of Betula (No. ha"^) 


StandDensityBemia 


doi:10.1371/journal.pone.0097192.t002 



the stand age of a Pinus koraiensis mixed forest, because the mean 
diameter of Pinus koraiensis in a stand naturally increases over time 
[16]. 

According to the recent CMIP5 (phase five of the Coupled 
Model Inter-comparison Project) climate change scenarios, the 
global temperature averaged in 2081-2100 is projected to likely 
exceed 2°C (RCP6.0 and RCP8.5 emission scenarios) above 
1850-1900; although unlikely to exceed 4°C in all other RCP 
emission scenarios than RCP8.5 [17]. The recent CMIP5 
scenarios also predicted that not only the surface air temperature 
but also the precipitation in 2071-2100 would be higher than 
those in 1986-2005 in Northeast China [18]. 

To examine the responses of composition of and structure of 
Pinus koraiensis mixed forests to spatial climatic changes, a field 
survey was conducted in the natural Pinus koraiensis mixed forests 
along a latitudinal gradient and an elevational gradient in 
Northeast China. The mathematical models were built for 
simulating the relationships of compositional and structural 
attributes of the Pinus koraiensis mixed forest with climatic and 
non-climatic factors, and two scenarios of the +2°C and +10% 
precipitation and +4°C and +10% precipitation were used to 
predict the possible changes of composition and structure in Pinus 
koraiensis mixed forests under the future global climatic change. 
The objectives of this study are to: (1) build mathematical models 
to quantify the responses of the composition and structure of Pinus 
koraiensis mixed forests to the climatic factors; (2) discover the key 
climatic factors which control the composition and structure of 
Pinus koraiensis mixed forest. Additionally, we hypothesized that the 
composition and structure of Pinus koraiensis mixed forests would 



significantiy be affected and Pinus koraiensis mixed forests would be 
overtaken by broadleaf-dominated forests in the region under the 
future global warming and increased precipitation. 

Materials and Methods 

Ethics Statement 

All necessary permits were obtained for the described field 
studies. This study was approved by State Key Laboratory of 
Vegetation and Environmental Change, Institute of Botany, the 
Chinese Academy of Sciences; Shengshan National Nature 
Reserve; Liangshui National Nature Reserve; Dongsheng Forestry 
Bureau; Changbai National Nature Reserve and Kuandian 
National Nature Reserve. 

Study Area 

This study was conducted in the Augusts of 2010 and 2012 in 
Northeast China, including ten sites: Kuandian, Dongsheng, 
Liangshui, Shengshan and Changbai from A to F (Fig. 1). The 
former four sites formed a latitudinal gradient while the latter six 
sites formed an elevational gradient. Each site consists of three 
plots (Table 1). All the 30 plots are natural Pinus koraiensis mixed 
forests without recent disturbance. AU the broadleaved trees are 
deciduous and most of the conifers are evergreen (except for Larix) 
in the study region. From June to August is the warmest quarter 
(summer) and January is the coldest month of the year in the study 
region. Rainy season begins at May or June and ends at 
September, and the summer (from June to August) precipitation 
occupies 60% of annual precipitation [19]. The 1 -km-resolution 
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bioclimate data of each plot were downloaded from WorldClim 
[20] (Table 2). 

Field Survey and Data Analyses 

In each plot, every living tree higher than 1 .4 m was identified 
to species and then measured for its circumference with a tape 
measure at its breast height (1.3 m). The measured circumference 
of each tree was converted into its diameter at breast height (DBH) 
and calculated into its basal area using the circle formulas. 

Fifteen attributes of Pinus koraiensis mixed forest were extracted 
from each plot. The fifteen attributes were; Mean DBH of Pinus 
koraiensis (MeanDBHp;n„s), Total Basal Area of Pinus koraiensis 
(BasalAreapiuus), Stand Density of Pinus koraiensis (StandDensityp;. 
nus)> Total Basal Area oi Acer and Fraxinus (BasalAreaAcer-Friixinus), 
Stand Density oiAcer and Fraxinus (StandDensityAcer-Friodnus), Total 
Basal Area of Quercus (BasalAreaQ^e„.^s), Stand Density of Quercus 
(StandDensityQuereus), Total Basal Area of Abies and Picea 
(BasalAreaAhies-Picea)) Stand Density of Abies and Picea (StandDen- 
sityAbie.s-Picea); Total Basal Area of Larix (BasalAreaLi^rix), Stand 
Density of Larix (StandDensityLarix), Total Basal Area of Ulmus 
(BasalAreauimus), Stand Density of Ulmus (StandDensityuimm), 
Total Basal Area of Betula (BasalAreaBetuia) and Stand Density of 
Betula (StancfDensityBetuia) (Table 2). 

One-way ANOVA and correlation and regression analyses were 
performed using SPSS 13.0. Standard errors within sites were 
detected by one-way ANOVA with descriptive options. Significant 
differences (p<0.05) between sites were detected using One-way 
ANOVA with post-Duncan's test. 

Stepwise linear regression analyses were performed in order to 
generate new mathematical models for the Pinus koraiensis mixed 
forest. The independent factors used in the regressions were linear 
and quadratic components of Max Temperature of Warmest 
Month (TempWarmestMonth), Min Temperature of Coldest 
Month (TempColdestMonth), Mean Temperature of Warmest 
Qriarter (TempWarmestQiiarter), Mean Temperature of Coldest 
Quarter (TempColdestQuarter), Precipitation of Wettest Month 
(PrecipWettestMonth), Precipitation of Driest Month (Precip- 
DriestMonth), Precipitation of Wettest Quarter (PrecipWettest- 
Quarter), Precipitation of Driest Quarter (PrecipDriestQuarter) 
and MeaiiDBHpinus. The 14 dependent variables were BasalAr- 
eapinus; StandDensitypinus, BasalAreaAcer-Fraxinus; StandDensityAcer- 
Friucinus, BasalAreaQj,^^^,, StandDensityQu^rcm, BasalAreaAbies-Pkca, 
StandDensityAhies-Picc-ii, BasalAreaLarix, StancfDensityLarix, BasalAr- 
eauimus, StandDensityuimus) BasalAreaBetuia and StandDensityB,,tuia 
(Table 2). Regression models with non-significant (p>0.05) term(s) 
were discarded. Among candidate models with all terms significant 
(p<0.05), the model with the lowest Small Sample Unbiased 
Akaike Information Criterion (AICc) value was selected [21]. 

Two-tailed partial correlation coefficients between each depen- 
dent variable and the independent variables excluded from its 
regression model were tested, controlling for the independent 
variable(s) included in its regression model. 

Running the Model 

The selected models were computed using the currently 
observed MeanDBHp;J,^s values under current climate, under 
the +2°C and +10% precipitation scenario and under the +4°C 
and +10% precipitation scenario. In order to examine the 
concordances between the simulated and observed values under 
current climate, scatter plots with simple error bars and regressions 
were generated using Sigmaplot 10.0. To represent the variations 
of stand densities and basal areas of Pinus koraiensis, other conifers 
and broadleaves under current climate under the +2°C and +10% 
precipitation scenario and under the +4°C and +10% precipitation 
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Table 5. Significant mathematical equations and regression coefficients (a, b, c, d) used to predict structure and composition of 
Pinus koraiensis mixed forests from linear and quadratic components of MeanDBHpinus, PrecipWettestMonth, TempWarmestQuater 
and TempWarmestMonth. 





Dependent 
Variable, Y 


Model Form 


Coefficients 












a (S.E.) 


b (S.E.) 


c (S.E.) d (S.E.) 




BasalAreapinus 


Y = a-l-b(MeanDBHpinus)^ 


14.740 (2.719)*** 


0.007 (0.002)** 




0.30** 


StandDensitypinus 


Y = a+b(PrecipWettestMonth)^+ 
c{MeanDBHpinus)+d{MeanDBHpinus)^ 


1232.331 (183.705)*** 


-0.003 (0.001)** 


-42.803 (9.983)*** 0.450 (0.136)** 


0.61*** 


BasalAreauimus 


Y = a+b(PrecipWettestMonth) 


-3.786 (1.222)** 


0.026(0.007)*** 




0.33*** 


StandDensityuimus 


Y = a+b(PrecipWettestMonth)^ 


-60.350(20.458)** 


0.003(0.001)*** 




0.45*** 


BasalAreaBetuia 


Y = a+b(PrecipWettestMonth) 


-6.445(2.200)** 


0.046(0.012)*** 




0.32*** 


StandDensitysetuia 


Y = a+b(PrecipWettestMonth) 


-129.583(23.427)*** 


0.874(0.129)*** 


_ _ 


0.61*** 


BasalAreaAcer-Fraxinus 


Y = a+b(TempWarmestQuater)+ 
c(MeanDBHpi„J^ 


-52.679(18.877)** 


3.018(1.057)** 


0.004(0.001)** 


0.29** 


StandDensityAcer-Fraxlnus 


Y = a+b{TempWarmestQuater) ^ 


-313.895(142.855)* 


1 .430(0.473)** 




0.22** 


BasalAreaQuercus 


Y = a+b(TempWarmestMonth) 


-161.720(57.358)** 


7.033(2.380)** 




0.21** 


StandDensityQuercus 


Y = a+b{TempWarmestQuater) 


-664.410(213.063)** 


41.112(12.278)** 




0.26** 


BasalAreaLarix 


Y = a+b(TempWarmestQuater) 


63.876(17.998)** 


-3.592(1.037)** 




0.28** 


StandDensityLarix 


Y = a-l-b{TempWarmestQuater) 


430.450 (110.937)*** 


-24.204 (6.393)*** 




0.32*** 


BasalAreaAbies-picea 


Y = a+b(TempWarmestMonth) 


141.160 (18.479)*** 


-5.679 (0.767)*** 




0.65*** 


StandDensityAbies-picea 


Y = a+b{TempWarmesti\/lonth) 


6330.509 (946.593)*** 


-254.657 (39.280)*** 




0.59*** 


Equations from stepwise regression analyses (n = 30). S.E, standard error; R^jj, adjusted multiple coefficient of determination. Significance level: (Non-significant) P> 
0.05. 



♦P<0.05. 

♦*P<0.01. 

***P<0.001. 

doi:l 0.1 371 /journal.pone.00971 92.t005 



scenario, vertical bar charts with error bars were generated using 
Sigmaplot 10.0. 

Results 

Structure and Composition of Pinus Koraiensis Forest 

Significant effects with latitude and/ or elevation were observed 
in most of the compositional and structural attributes of Pinus 
koraiensis mixed forest (P<0.05, Table 3). The StandDensitypinus 
significantiy increased while BasalAreauimusj StandDensityuimu,,, 
BasalAreaBetuia, StandDensityBetuia, BasalAreaAcer-Fraxinm and 
StandDensityAcer-Fraxinus significantly declined with increasing 
latitude (from the southernmost site Kuandian to the northern- 
most site Shengshan) (P<0.05, Table 3). BasalAreaLa^jx, Stand- 
DensityLarix; BasalAreaAbies-picea and StandDensityAbies-Picea were 
significandy enhanced with increasing elevation (from Changbai A 
to Changbai F) (P<0.05, Table 3). 

Responses of Composition and Structure in Pinus 
Koraiensis Forests to Climatic and Non-Climatic Factors 

A few excluded independent variables showed significant partial 
correlations with the dependent variable (P<0.05, Table 4). The 
significant negative relationships of the BasalAreapinus with the 
climatic factors of PrecipWettestMonth and PrecipWettestQuar- 
ter, of the BasalAreaBetuia and StandDensityBetuia with the climatic 
factors of PrecipDriestMonth and PrecipDriestQjaarter, of the 
BasalAreaBetuia with the climatic factor of TempColdestMonth, 
and of the StandDensityBetuia with the non-climatic factor of 
MeanDBHpinus were observed (P<0.05, Table 4). 



Significant regression models were observed for all the 14 
dependent variables (P<0.05, Table 5). MeanDBHp;nu., had 
positive relationships with BasalAreapinus and BasalAreaAcer- 
Fraxinus, aiid had negative relationship with StandDensitypinus 
when the MeanDBHjJinus was below 47.6 cm. PrecipWettest- 
Month suppressed StandDensitypinus, while promoted BasalAr- 
eauimus, StandDensityuimus, BasalAreaBetuia and StandDensityBe- 
tuia- TempWarmestQuarter enhanced BasalAreaAeer-Fraxinus, 
StandDensityAcer- Fraxinm and StandDensityquercus, but suppressed 
BasalAreaLarix and StandDensityLarix- TempWarmestMonth pro- 
moted BasalAreaQuercus, but suppressed BasalAreaAbie.s-Kcea and 
StandDensityAbies-Picea (Table 5). 

Model Outputs 

There were better concordances between observed and 
simulated values in stand densities than in basal areas. More 
significant linear regression relationships between observed and 
simulated values with higher adjusted R squares and closer-to-1 
slopes were found in StandDensitypinus (Fig. 2a), StandDensity- 
Quercus (Fig. 2c), StandDcusityuimus (Fig. 2d) and StandDensityBe- 
tuia (Fig- 2e), compared with those in BasalAreapinus (Fig- 2h), 
BasalAreaQuercus (Fig. 2j), BasalAreauimus (Fig- 2k) and BasalAr- 
eaBetuia (Fig- 21)- However, better linear regression relationships 
between simulated and observed values were found in BasalAr- 
eaAcer-Fraxinus (Fig- 21) and BasalAreaAbies-Picea (Fig- 2m), compared 
with those in StandDensityAcer-Fraxinus (Fig- 2b) and StandDensi- 
tyAbies-Picea (Fig. 21). The linear regression relationships between 
simulated and observed values in StandDensityLarbc (Fig- 2g) and 
BasalAreaLarix (Fig- 2n) were not significant. 
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Figure 2. Linear regression relationships between the simulated (n = 3, mean with standard error) and observed (mean of each site) 
values of StandDensitypr„us (No. ha \- a), StandDensityAcer-Fraxfnus (No. ha'^; b), StandDensityQ„ercu5 (No. ha ^ c), StandDensi- 
tyuimus (No.ha"'; d), StandDensitysetuia (No.ha \' e), StandDensityAbfes-pkea (No. ha \- f), StandDensityLarix (No. ha ^ g), 
BasalAreapinus (m^ ha \- h), BasalAreaA„,_F„xinus (m^ ha \- i), 8asalAreaQ„ercus (m^ ha"^; j), BasalAreauimus (m^ ha k), 
BasalAreaBetuia (m^ ha I), BasalAreaAbies-Pkea (n™^ ha~'; m) and BasalAreaLarix (m^ ha~'; n) in Pinus koraiensis mixed forests. 
Significance level: '~'^(Non-significant)P>0.05, *P<0.05, **P<0.01, ***P<0.001. 
doi:1 0.1 371/journal.pone.00971 92.g002 
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Figure 3. Stand density of Pinus koraiensis (a), total basal area of Pinus Icoraiensis (b), stand density of other conifers (c), total basal 
area of other conifers (d), stand density of broadleaves (e), total basal area of broadleaves (f) of Pinus icoraiensis and other conifers 
[Abies and Picea and Larij^l and broadleaf trees (broadleaves) In the ten sites (from left to right: Kuandlan, Dongsheng, Llangshul, 
Shengshan and Changbai from A to F) observed under current condition and modeled under current climate and under +2°C and 
+10% precipitation scenario and under +4°C and +10% precipitation scenario (using the currently observed IVleanDBHpi„us values). 
doi:1 0.1 371/journal.pone.0097192.g003 



Stand densities and total basal areas of conifers were predicted 
to decrease wliUe stand density and total basal area of broadleaved 
trees were predicted to increase under the +2°C and +10% 
precipitation and +4°C and +10% precipitation scenarios. Stand 
density of Pinus koraimsis was predicted to notably decrease in 
Kuandian and slighdy decrease in other sites (Fig. 3a). Total basal 
area of Hnus koraiensis was predicted to keep constant (Fig. 3b). 
Stand density and total basal area of other conifers {Abies and Picea 
and Larix) were predicted to be zero (Fig. 3c and d). Broadleaved 
trees would increase dramatically (Fig. 3e and f). In all the ten sites, 
the densities and total basal areas of broadleaved trees would be 
higher than those of Pinus koraiensis under the +2°C and +10% 
precipitation and +4"C and +10% precipitation scenarios (Fig. 3). 

Under each climate scenario, the stand density of Pinus koraiensis 
(StandDensitypi,„,s) firstly increased and then decreased with 
increasing latitude, whUe the stand density of broadleaved trees 
decreased with increasing latitude and altitude (Fig. 3). 

Discussion 

Betula and Ulmus serve as early successional species (pioneers), 
while Aeer and Fraxinus and Pinus koraiensis serve as late successional 
species during the succession oi Pinus koraiensis mixed forest [1,22]. 
In our results, MeanDBHpjjms (which represents the stand age of a 
Pinus koraiensis forest) [16] was observed to have positive 
relationships with BasalAreapinus and BasalAreaAcer-Frajdnus) sug- 
gesting that Pinus koraiensis and Acer and Fraxinus are late 
successional species and naturally increase in their total basal 
areas o\'cr tim(;. However, negative relationship <;xisti'd between 
MeanDBHpi„us and StandDensitypmus when MeanDBHpjjms was 
below 47.6 cm (Table 5), which indicated the self-thinning 
phenomenon oi Pinus [23]. 

Heavier rain usually results in higher soil wetness [24]. Soil 
wetness greatly reduces soil cohesion and increases the possibility 
of tree uprooting [25]. Uprooting of trees results in exposed 
mineral soil (bare soil) [26,27] . The exposed mineral soil favors the 
regeneration of pioneer trees such as Betula [28] and Ulmus [29]. In 
our results, PrecipWettestMonth was observed to have negative 
effect on StandDensitypi„us, but positive effects on BasaIAreauijj,u„ 

StandDensityuimm, BasalAreaBetuia and StandDensityBetuia 
(Table ,5). In addition, better concordances between simulated 
and observed values were found in StandDensitypinu,, StandDen- 
sityuim and StandDensityBetuia, compared with those in BasalAr- 
eapinus, BasaLAreauim and BasalAreaBetuia (Fig- 2). Our results 
indicate that higher monthly precipitations can cause higher soil 
wetnesses, more uprooting events and more exposed mineral soil, 
and therefore negatively modulate the StandDensitypinus and 
positively modulate the stand densities of Betula and Ulmus. 

CUmate warming has positive and negative elfects on plants. 
Prolonged growing seasons and increased growing-season temper- 
atures have been found to enhance the growth of some plants, 
especially the growth of many deciduous broadleaved trees 
[30,31]. But during midday of summer, high temperature can 
inhibit the photosynthesis and growth of plants by provoking leaf- 
to-arr vapor pressure deficit (VPD) [32,33,34]. Generally, broad- 
leaf trees are less sensitive than conifers trees [35] while Pinus trees 
are less sensitive than other conifer trees {Picea and Larix) [36,37] to 
the leaf-to-air VPDs. In our result, no effects of temperature on 



Pinus koraiensis but the negative effects of summer temperatures on 
other conifers {Abies and/ or Picea and Larix) and the positive effects 
of summer temperatures on some broadleaved trees {Acer and/ or 
Fraxinus and Quercus) were observed (Table 5). Our results indicate 
that (1) hotter summers have hotter middays and thus more severe 
leaf-to-air VPDs, which are more beneficial to the growth of 
broadleaf trees compared with conifer trees, and more beneficial 
to the growth of Pinus koraiensis trees compared with other conifer 
trees {Abies and/ or Picea and Larix); (2) increased summer (growing- 
season) temperatures may enhance the growths of many broad- 
leaved trees {Acer and/ or Fraxinus and Quercus). 

A wide-ranged pine species is usually composed of many 
genetically-different local populations. Each local population has 
evolved to adapt to its local climate, and has its unique response 
curve to climate change [11,38]. In our results, no effect of 
temperature on the density and total basal area of Pinus koraiensis 
were observed (Table 5), and the reason may possibly is due to the 
local adaptations of many geneticaUy-drfiFerent Pinus koraiensis 
populations. We suggest that the Pinus koraiensis in the region may 
not necessarily maintain itself under the +2"C and +10% 
precipitation and +4°C and +10% precipitation scenarios if 
proper migrations of local populations do not occur in time. 

Studies using a forest gap model alone usually predicted the 
extinction of I^nus koraiensis within 150 years in current L^nus 
koraiensis mixed forest under climate warming scenarios [7,39], 
while other study linking a forest gap model with a landscape 
model predicted that the Pinus koraiensis could persist for at least 
300 years under the warming climate [7]. Our result agreed with 
the latter by predicting that F^nus koraiensis would persist, while the 
broadleaf trees would o\'ertake Pinus koraiensis to become dominant 
species under either of the two warming scenarios (Fig. 3). 

Conclusions 

These results show that some mathematical models were 
successfully built to quantify the responses of the composition 
and structure of Pinus koraiensis mbced forests to the climatic factors. 
TempWarmestMonth, TempWarmestQuarter and PrecipWettest- 
Month strongly affected the stand densities and total basal areas in 
Pinus koraiensis mixed forest and therefore they are suggested to be 
key climatic factors which shape the composition and structure of 
Pinus koraiensis mixed forest. The composition and structure oi Pinus 
koraiensis mixed forests were significantiy affected, and although 
Pinus koraiensis would persist, the current forests dominated by Pinus 
koraiensis in the region would all shift and become broadleaf- 
dominated forests due to the dramatic increase of broadleaf trees 
under the future global warming and increased precipitations. 
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